Weighted evolving networks: coupling topology and weights dynamics 
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We propose a model for the growth of weighted networks that couples the establishment of new 
edges and vertices and the weights' dynamical evolution. The model is based on a simple weight- 
driven dynamics and generates networks exhibiting the statistical properties observed in several 
real-world systems. In particular, the model yields a non-trivial time evolution of vertices' properties 
and scale-free behavior for the weight, strength and degree distributions. 
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The last years have witnessed a hectic activity devoted 
to the characterization and understanding of networked 
structures as diverse as ecological and biological systems 
or the Internet and the WWW Q, H H H- These net- 
works generally exhibit complex topological properties 
such as the small-world phenomenon [5| and scale-free 
behavior Q . The need for explaining these complex topo- 
logical features has led to a new modeling framework that 
focuses on the dynamical evolution and growth of net- 
works 0. In this perspective, a wide array of models 
aimed at capturing various properties of real networks 
have been formulated Q. These models, however, do 
generally consider only the topological structure and do 
not take into account the interaction strength -the weight 
of the link- that characterizes real networks. For exam- 
ple, the diversity of the predator-prey interaction and 
of metabolic reactions is considered a critical ingredient 
of ecosystems 0, 13 and metabolic networks , respec- 
tively; in social systems, the weight of interactions is very 
important in the characterization of the corresponding 
networks 0] • Similarly, the Internet traffic J3 or the 
number of passengers in the airline network [zllll 
are crucial quantities in the study of these systems. In- 
terestingly, recent studies j3, liii LLJ] have shown that 
weighted networks exhibit additional complex properties 
such as broad distributions and non-trivial correlations 
of weights that do not find an explanation just in terms 
of the underlying topological structure. 

In this Letter, we define a simple model for weighted 
networks that considers the basic evolution of the system 
as driven by the weight properties of edges and vertices. 
In addition, differently from early models proposed in the 
past |15lll6| |. we allow the dynamical evolution of weights 
during the growth of the system. This mimics the evo- 
lution and reinforcements of interactions in natural and 
infrastructure networks. The generated networks display 
power-law behavior for the weight, degree and strength 
distributions with non-trivial exponents depending on 
the unique parameter defining the model's dynamics. 
These results suggest that the inclusion of weights in 
networks modeling naturally explains the diversity of 



scale- free behavior empirically observed in real networked 
structures. Strikingly, the weight-driven growth recovers 
an effective preferential attachment for the topological 
properties, providing a microscopic explanation for the 
ubiquitous presence of this mechanism. 

Weighted networks are usually described by a matrix 
Wij specifying the weight on the edge connecting the ver- 
tices i and j , with i,j = 1 , . . . , N where N is the size of the 
network (wij = if the nodes i and j are not connected). 
In the following we will consider only the case of symmet- 



ric weights Wi 



Prototypical examples of weighted 



networks can be found in the world-wide airport network 
and the scientific collaboration network 
(SCN) |l7l lSj . In the airport network each given weight 
Wij is the number of available seats on direct flights con- 
nections between the airports i and j. In the SCN the 
nodes are identified with authors and the weight depends 
on the number of co-authored papers 0, 0. A first 
topological characterization of networks is obtained with 
the vertices degree k, number of connected neighbors, 
and the associated probability distribution P(k). Com- 
plex networks often exhibits a power-law degree distribu- 
tion P(k) ~ fc -7 with 2 < 7 < 3 and measurements of 
weighted networks confirm this evidence 

00 Sim- 
ilarly, a first characterization of weights is obtained by 
the distribution P(w) that any given edge has weight w. 
In complex networks such as the WAN and SCN this 
distribution is heavy tailed and spans several orders of 
magnitude 0,0. Along with the degree of a node, a 
very significative measure of the network properties in 
terms of the actual weights is obtained by looking at the 
vertex strength Sj defined as [T3 . Il5| 



jev(i) 



'ij-, 



(1) 



where the sum runs over the set V[i) of neighbors of i. 
The strength of a node integrates the information about 
its connectivity and the weights of its links, and can be 
considered as the natural generalization of the connec- 
tivity. For instance, in the case of the WAN the strength 
provides the actual traffic going through the vertex i and 
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FIG. 1: Illustration of the construction rule. A new node 
n connects to a node i with probability proportional to 
Si I y\ Sj . The weight of the new edge is wo and the total 
weight on the existing edges connected to i is modified by an 
amount equal to 5. 

is an obvious measure of the size and importance of each 
airport. This quantity obviously increases with the ver- 
tex degree ki and usually displays the power-law behavior 
s ~ k@ , with the exponent (3 varying as a function of the 
specific network. Finally, we note that in most cases, the 
empirical distribution of s is heavy-tailed analo- 
gously to the power-law decay of the degree distribution. 

Previous approaches to the modeling of weighted net- 
works focused on growing topologies where weights were 
assigned statically, i.e. once for ever, with different rules 
related to the underlying topology [ijj |lj| . These mech- 
anisms, however, overlook the dynamical evolution of 
weights according to the topological variations. We can 
illustrate this point in the case of the airline network. If 
a new airline connection is created between two airports 
it will generally provoke a modification of the existing 
traffic of both airports. In general, it will increase the 
traffic activity depending on the specific nature of the 
network and on the local dynamics. In the following we 
will implement a model that takes into account the cou- 
pled evolution in time of topology and weights. 

The model dynamics starts from an initial seed of N 
vertices connected by links with assigned weight w - At 
each time step, a new vertex n is added with m edges that 
are randomly attached to a previously existing vertex i 
according to the probability distribution 

n„_ = (2) 

This rule relaxes the usual degree preferential attach- 
ment, focusing on a strength driven attachment in which 
new vertices connect more likely to vertices handling 
larger weights and which are more central in terms of 
the strength of interactions. The weight of each new 
edge is fixed to a value Wq. Moreover, the presence of the 
new edge (n, i) will introduce variations of the existing 
weights across the network. In particular, we consider 
the local rearrangements of weights between i and its 
neighbors j e V(i) according to the simple rule 

Wij — > Wij + Awi 3 , (3) 



where 

A Wij = 5^i. (4) 

Si 

This rule considers that the establishment of a new edge 
of weight wq with the vertex i induces a total increase 
of traffic 8 that is proportionally distributed among 
the edges departing from the vertex according to their 
weights (see Fig. HJ, yielding s, — > + 8 + wq. At this 
stage, it is worth remarking that while we will focus on 
the simplest model with 8 = const, different choices of 
Amy with heterogeneous Si or depending on the spe- 
cific prop erties of each vertex (toy, ki, s,) can be consid- 
ered [lj(. Finally, after the weights have been updated 
the growth process is iterated by introducing a new ver- 
tex with the corresponding re-arrangement of weights. 

The weight driven attachment (Eq. (J2J) appears to be a 
plausible mechanism in many technological networks. In 
the Internet new routers connect to more central routers 
in terms of bandwidth and traffic handling capabilities 
and in the airport networks new connections are gener- 
ally established to airports with a large passenger traffic. 
Even in the SCN this mechanism might play a role since 
an author with more co-authored papers is more visible 
and open to further collaborations. Also the weights dy- 
namics can be traced back to simple dynamical processes. 
Indeed, the quantity u> sets the scale of the weights and 
without loss of generality we can use the rescaled quan- 
tities / wo, Si I 'wq and 8/wq. For the sake of simplicity 
we can set wq = 1 and the model depends only on the 
dimensionless parameter 8, that is the fraction of weight 
which is 'induced' by the new edge onto the others. Ac- 
cording to the value of 8, different scenarios are possible. 
If the induced weight is 8 w wo — 1 we mimic situations 
in which an appreciable fraction of traffic generated by 
the new connection will be dispatched in the already ex- 
isting connections. This is plausible in the airport net- 
works where the transit traffic is rather relevant in hubs. 
In the case of 8 < I we face situations such as the SCN 
where it is reasonable to consider that the birth of a new 
collaboration (co-authorship) is not triggering a more in- 
tense activity on previous collaborations. Finally, 8 > 1 
is an extreme case in which a new edge generates a sort of 
multiplicative effect that is bursting the weight or traffic 
on neighbors. 

In the model, the time is measured with respect to the 
number of vertices added to the graph, resulting in the 
definition t = N — Nq, and the natural time scale of the 
model dynamics is the network size N. The network's 
evolution can be inspected analytically by studying the 
time evolution of the average value of Sj (t) and ki (t) of 
the i-th. vertex at time t, and by relying on the contin- 
uous approximation that treats k, s and the time t as 
continuous variables 0, Q ■ When a new edge is added 
to the network, the strength Sj of vertex i can increase 
either if the edge connects directly to i or to one of its 
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FIG. 2: (A) Evolution of the strength of vertices during the 
growth of the network, for various values of 5; The thick lines 
are the predicted power laws t a , a = (l + 2<5)/(2 + 2<5) (m = 2, 
N = 10 4 ). (B) Strength s t versus ki for various values of 8; 
the straight lines are the predictions Sj = (1+ 2S)ki. 



neighbor. The evolution equation for s, is thus given by 

dSi Sj ,^ „^ x > Sn 

~dt 



28 



1 «i(t) 



si 



28 + 2 t 



(5) 



where the last expression is recovered by noticing that 
each added link increases the total strength by an amount 
equal to 2 + 26 thus implying that Si(t) w 2m(l + <5)t. 
This equation can be readily integrated with the initial 
condition Si(t = i) = m, yielding 



»(*) = 



2 A j- 1 
2 rV -|- 2 



(6) 



Similarly, the rate equation for the degree evolution reads 
dki Si(t) s l (t) 



dt 



E,*i(*) 2(1 + «5)t : 



(7) 



that implies fci(t) = Si(t)/(26 + 1). The proportionality 
relation s ~ fe, that implies /? = 1, is particularly relevant 
since it states that the weight-driven dynamics generates 
in Eq. an effective degree preferential attachment that 
is parameter independent. This highlights an alternative 
microscopic mechanism accounting for the presence of the 
preferential attachment dynamics in growing networks. 

The behavior of the various statistical distributions can 
be computed by noticing that the time at which the node 
i enters the network is uniformly distributed in [0,t] and 



P(M) 



l 



t + N J 



6{k - ki(t))di, 



(8) 



where S(x) is the Dirac delta function. Solving the above 
integral and considering the infinite size limit t — > oo 
yields P(k) ~ fc~ 7 with 

7=^ (9) 
7 26+ 1 () 
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FIG. 3: (A) Probability distribution P(s). Data fitting points 
(filled circles) are consistent with the analytic prediction of 
power-law behavior s -7 with 7 = (3 + 4(5) /(l + 28) (solid 
line) as shown in the inset. (B) Probability distribution of 
the weights P{w) ~ w~ a . In the inset we report the value 
of a obtained by data fitting (filled circles) and the analytic 
expression a — 2 + 1/5 (solid line). The data are averaged 
over 200 networks of size N = 10 4 . 



It is straightforward to obtain the same behavior P(s) ~ 
s 1 for the strength distribution since s and k are pro- 
portional. This result shows that the obtained graph is 
a scale- free network described by an exponent 7 G [2,3] 
that depends on the value of the parameter 5. In par- 
ticular, when the addition of a new edge doesn't affect 
the existing weights (6 = 0), the model is topologically 
equivalent to the Barabasi- Albert model |fj and the value 
7 = 3 is recovered. For larger values of 6 the distribu- 
tion is progressively broader with 7-^2 when 6 — > 00. 
This indicates that the weight-driven growth generates 
scale-free networks with exponents varying in the range 
of values usually observed in the empirical analysis of net- 
worked structures Q, 0, • Noticeably the exponents are 
non-universal and depend only on the parameter 6 gov- 
erning the microscopic dynamics of weights. The model 
therefore proposes a general mechanism for the occur- 
rence of varying power-law behaviors without resorting 
on more complicate topological rules and variations of 
the basic preferential attachment mechanism. 

In order to check the analytical predictions we per- 
formed numerical simulations of networks generated by 
using the present model with different value of 6, mini- 
mum degree m and varying size N. In Fig. we show the 
behavior of the vertices' strength versus time for different 
value of 6, recovering the behavior predicted analytically. 
We also report the average strength Si of vertices with 
degree ki and confirm that /? = 1 as well as the value 
of the pre- factor 26 + 1. Finally, in Fig.^K we plot the 
probability distribution P(s) that exhibits a power-law 
behavior in excellent agreement with the theoretical pre- 
dictions. 

Similarly to the previous quantities, it is possible to 
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FIG. 4: Time evolution of Wij during the growth of the net- 
work (m = 2 and N = 10 4 ) for different values of <5. The func- 
tional behavior is consistent with the predicted power laws t b , 
b = 5/(1 + S), shown by thick lines. 



properties, a natural generalization of the model consists 
in considering 8 as a function of the vertices degree or 
strength. For example, a non-linear generalization of eq. 
P)l can yield non-linear correlations between strengths 
and degrees (/3 ^ 1) 0]. Similarly, more complicated 
variations of the microscopic rules may be implemented 
to mimic in a detailed fashion particular networked sys- 
tems. In this perspective the present model appears as a 
general starting point for the realistic modeling of com- 
plex weighted networks. 
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by the EC-Fet project COSIN IST-2001-33555. 



obtain analytical expressions for the evolution of weights 
and the relative statistical distribution. The weight Wij 
increases by the addition of a new link either on i or on 
j and the corresponding rate equation can be written as 



dwij 



E 



3 S i 



8^- 

Si 



j s j s j 



(10) 



This equation can be solved with the initial condition 
Wij(tij) — 1, where Uj =max(i,j) is the time at which 
the edge is established, yielding Wij(t) — (t/t i:j ■) s / t - s+1 ^ 
From this expression it is possible to calculate similarly 
to P(k) the probability distribution P(w) which is in this 
case also a power-law P(w) ~ w~ a where 



(11) 



The exponent a has large variations as a function of the 
parameter 8 and P{w) moves from a delta function for 
8 = to a very slow decaying power-law with a = 2 
if 8 — > oo. This feature clearly shows that the weight 
distribution is extremely sensitive to changes in the mi- 
croscopic dynamics ruling the network's growth. These 
predictions can be confirmed by numerical simulations 
and the results are shown in Figs.|3}3 and0] 

The model we have presented is possibly the simplest 
one in the class of weight-driven growing networks. A 
novel feature in the model is the weight dynamical evo- 
lution occurring when new vertices and edges are intro- 
duced in the system. This simple mechanism produces 
a wide variety of complex and scale-free behavior de- 
pending on the physical parameter 8 that controls the 
local microscopic dynamics. While a constant parameter 
8 is enough to produce a wealth of interesting network 
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